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Results of theoretical calculations of ionization rates of Rb and Na Rydberg atoms by blackbody 
radiation (BBR) are presented. Calculations have been performed for nS, nP and nD states of 
Na and Rb, which are commonly used in a variety of experiments, at principal quantum numbers 
n=8-65 and at three ambient temperatures of 77, 300 and 600 K. A peculiarity of our calculations 
is that we take into account the contributions of BBR-induced redistribution of population between 
Rydberg states prior to photoionization and field ionization by extraction electric field pulses. The 
obtained results show that these phenomena affect both the magnitude of measured ionization rates 
and shapes of their dependences on n. The calculated ionization rates are compared with the results 
of our earlier measurements of BBR-induced ionization rates of Na nS and nD Rydberg states with 
n=8-20 at 300 K. A good agreement for all states except nS with n > 15 is observed. We also present 
the useful analytical formulae for quick estimation of BBR ionization rates of Rydberg atoms. 

PACS numbers: 32.80.Fb, 32.80.Rm, 32.70.Cs 



I. INTRODUCTION. 

Blackbody radiation (BBR) is known to strongly af- 
fect the populations of atoms in highly excited Ryd- 
berg states It has also been shown that at the am- 
bient temperature of 300 K BBR can photoionize Ry- 
dberg atoms with n '--^ 20 at astonishingly high rates 
(~ 103 s-i) [1, i, II. Strong effect of BBR on Rydberg 
atoms is related to large matrix elements of bound-bound 
and bound-free transitions between Rydberg states in the 
microwave and far infrared spectral range @ . 

Interaction of Rydberg atoms with BBR has been stud- 
ied earlier in various contexts. Farley and Wing [6i] cal- 
culated the dynamic Stark shifts and depopulation rates 
of Rydberg levels of alkali atoms with n < 30 at 300 K. 
The temperature dependence of BBR-induced transitions 
rates from the 19S state of sodium was calculated and 
measured by Spencer et al [3] . Galvez et al 0, [1] have 
studied BBR-induced cascade transitions from the ini- 
tially populated n—2A-29 states of Na, both theoretically 
and experimentally. 

Although interaction of Rydberg atoms with black- 
body radiation has been studied for years, both theoreti- 
cally and experimentally, only a few works were devoted 
to BBR-induced ionization of Rydberg atoms itself. The 
temperature dependence of the BBR-induced ionization 
rate of Na 17D state was numerically calculated using 
the quantum defect method in a Coulomb approxima- 



tion and experimentally measured by Spencer et al 
A simple scaling law for BBR ionization rates was also 
introduced in that work. More recently, the interest to 
BBR-induced ionization of Rydberg atoms has been re- 
lated to the spontaneous formation of ultracold plasma 
in dense samples of cold Rydberg atoms [l3| , and to 
the prospects of its use as a convenient reference signal in 
absolute measurements of coUisional ionization rates [11 1 . 
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Even nowadays, however, the studies of spontaneous evo- 
lution of ultracold Rydberg atoms to a plasma caused 
by BBR I, [13] use the simple estimates of BBR ioniza- 
tion rates taken from the well-known work The nu- 
merical calculations of direct BBR photoionization rates 
were performed by Lehman ^T5| for H and Na Rydberg 
states with n= 10-40 using a Herman-Skillman potential 
with an additional term to account for core polarization. 
However, the theoretical data for Rb atoms, which are 
widely investigated in experiments on ultracold plasma, 
are lacking. No systematic experimental studies of the 
n-dependences of BBR ionization rates of alkali-metal 
Rydberg atoms in a wide range of principal quantum 
numbers are known. 

A major problem in the interpretation of measured 
ionization rates of Rydberg atoms is related to the 
fact that populations of Rydberg states are affected by 
BBR-induced processes to an unknown (or, at least, 
not straightforwardly predictable) extent, which de- 
pends on the combination of specific experimental condi- 
tions (principal quantum numbers n, ambient tempera- 
ture, duration of measurements, extracting electric field 
strengths, etc.). Therefore, a more detailed study of 
BBR-induced ionization of Rydberg atoms under typical 
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experimental conditions is required, especially for higher 
Rydberg states that are often explored in the experiments 
with cold atoms. 

In this article we present the results of numerical cal- 
culations of BBR ionization rates of Rb and Na Ryd- 
berg atoms for the most commonly used nS, nP and nD 
states with n=8-65 at the ambient temperatures of 77, 
300 and 600 K. In these calculations we take into account 
two phenomena that may affect the observed ionization 
rates: the time-dependent BBR-induced population re- 
distribution between Rydberg states prior to photoion- 
ization and the selective field ionization (SFI) of high- 
lying Rydberg states by the electric field applied for ex- 
traction of ions from the excitation zone. Although both 
phenomena were mentioned in some of the earlier stud- 
ies 01 , their effects have not yet been studied in sufficient 
detail neither theoretically, nor experimentally. In partic- 
ular, in the measurements of direct BBR photoionization 
rate of the sodium 17D state 0] with effective lifetime 
4 fis at the ambient temperature of 300 K, the effect of 
population redistribution was diminished due to the short 
measurement time of 500 ns and weak electric field used 
for extraction of ions. However, for Rydberg states with 
n 8 and effective lifetimes less than 500 ns even this 
very short time interval is insufficient to avoid mixing 
with neighboring states, so that mixing processes must 
be necessarily accounted in calculations. In the present 
study we address to both phenomena and provide the 
useful analytical formulae, which can be applied to any 
Rydberg atom for estimates of the direct BBR-induced 
ionization rates and of the contribution of the SFI. 

The simplest (but often insufficient) way of considering 
BBR-induced ionization after the excitation of an atom 
^ to a given nL Rydberg state is to limit the problem to 
the direct photoionization of the initial nL state in one 
step by absorption of BBR photons: 

A{nL) + nujBBR.^ A+ +e-, (1) 

where fvjjBBR is the energy of absorbed BBR photon, 
is the atomic ion, and e~ is the free electron emitted 
in the ionization. In the reality, however, ionization of 
Rydberg atoms exposed to BBR is a complex process, 
in which the following main components can be identi- 
fied [see Fig. 1(a)]: (i) direct photoionization of atoms 
from the initial Rydberg state via absorption of BBR 
photons, (ii) field ionization by extraction electric field 
pulses of high Rydberg states, which are populated from 
the initial Rydberg state by absorption of BBR photons, 
(iii) direct BBR-induced photoionization of atoms in the 
neighboring Rydberg states, which are populated due to 
absorption and emission of BBR photons prior to pho- 
toionization, and (iv) field ionization of other high-lying 
states, which are populated via population redistribution 
involving two or more steps of BBR photon absorption 
and/or emission events. Our calculations show that all 
these processes can contribute to the total ionization rate 
to a comparable extent, and, therefore, none of them 
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FIG. 1: (a) Schematic illustration of BBR induced and field 
ionization processes occurred after excitation of the initial 
Na(16S) state, including redistribution of population over 
other n'L' Rydberg states due to spontaneous and BBR in- 
duced transitions from the initial state, (b) Calculated spon- 
taneous and BBR induced transition rates from the initial 
16S state to other n'P states, (c) Calculated spontaneous 
and BBR induced transition rates from the initial 16D state 
to other n'P states. 

can be disregarded. In what follows we will consider the 
above processes separately and calculate the total BBR 
ionization rates, both analytically and numerically. 

II. CALCULATION OF BBR IONIZATION 
RATES. 

Ionization mechanisms of Rydberg atoms exposed to 
BBR are illustrated in Fig. 1. The total BRR- induced 
ionization rate can be written as consisting of four sepa- 
rable contributions: 

W'b%r = Wbbr + Wsfi + M/BBfl + WIsTi ■ (2) 

The first contribution, Wbbr, is the direct BBR pho- 
toionization rate of the initially excited nL state, which 
will be discussed in Section II A. The second term, 
WsFi, is the rate of the selective field ionization (SFI) 
of high n"L' Rydberg states, which are populated from 
the initial Rydberg state nL via absorption of BBR pho- 
tons. This field ionization will be discussed in Sec- 
tion II C, while redistribution of population between Ry- 
dberg states will be described in Section II B. The third 
term, Wbjbr, is the total rate of BBR-induced photoion- 
ization of neighboring n' L' Rydberg states, which are 
populated via spontaneous and BBR-induced transitions 
from the initial state. The last term, Wg]?f , is the rate 
of SFI of high-lying Rydberg n" L' states that are popu- 
lated in a two-step process via absorption of BBR pho- 
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tons by atoms in n'L' states (note, that here, in contrast 
to WsFi, we consider lower n'L', states which cannot be 
directly field ionized). These latter two ionization rates, 
which are related to population redistribution between 
Rydberg states, will be considered in Section II D . The 
atomic units are used below, unless specified otherwise. 



rate depends on parameters n, L, and T. Such ex- 
pressions can be obtained using the analytical formu- 
lae for bound-bound and bound-free matrix elements 
deduced by Goreslavsky, Delone and Krainov (GDK) 
(l^ in the quasiclassical approximation. For the direct 
BBR-induced photoionization of an nL Rydberg state the 
cross-section is given by: 



A. Direct BBR photoionization. 

Direct BBR-induced photoionization rate Wbbr of a 
jivcn nL state is calculated from the general formula [2| : 



a^{nL ^ E,L±1) 



LUL" 



(6) 



Wbbr = c p^duj, 



(3) 



where c is the speed of light, ujnL = 1/ (^^rtgyjj is the 

photoionization threshold frequency for the nL Rydberg 
state with the effective principal quantum number n^f / = 
(n — 6l), where Sl is a quantum defect, and <7uj is the 
photoionization cross-section at the frequency lu. The 
volume density p^^ of BBR photons at the temperature 
T is given by the Plank distribution: 



7r2c3 [e'^/ikT) _ 1] 



(4) 



where kT is thermal energy in atomic units. For 
isotropic and non-polarized thermal radiation field, the 
value of (T^ is determined by the radial matrix elements 
R {nL ^ E, L ±1) oi dipole transitions from discrete nL 
Rydberg states to the continuum states with L±l and 
photoelectron energy E: 



^c{2L+l) ^^^^ 



LraaxR'inL^ E,L±l), (5) 



where Lmax is the largest of L and L'. 

The main problem in the calculation of Wbbr for 
an arbitrary Rydberg state is thus to find the values 
of R{nL^E,Lztl) and their frequency dependence. 
In order to achieve high accuracy of the matrix ele- 
ments, numerical calculations should be used. In this 
work we used the semi-classical formulae derived by Dy- 
achkov and Pankratov Il3l . In comparison with other 
semi-classical methods [l4L [lH. this method is advanta- 
geous as it gives orthogonal and normalized continuum 
wavefunctions, which allow for the calculation of pho- 
toionization cross-sections with high accuracy. We have 
verified that photoionization cross-sections of the lower 
sodium S states calculated using the approach of [l^l 
are in good agreement with the so phi sticated quantum- 
mechanical calculations by Aymar |16j . 

Approximate analytical expressions for Wbbr would 
also be useful, since they illustrate how the ionization 



where K,^ (x) is the modified Bessel function of the sec- 
ond kind. This formula was initially derived to describe 
the photoionization of hydrogen atom. It was assumed 
that the formula can be extended to alkali atoms sim- 
ply by replacing n by n^ff = {n — Sl), where Sl is the 
quantum defect of the Rydberg state. In the reality, how- 
ever, its accuracy in absolute values is acceptable only for 
truly hydrogen- like states with small quantum defects. A 
disadvantage of the GDK model is that it disregards non- 
hydrogenic phase factors in the overlap integrals of dipole 
matrix elements. 

The main contribution to Wbbr in Eq. ^ is due to 
BBR frequencies near the ionization threshold frequency 
i^nL , because the bound- free dipole moments rapidly de- 
crease with increasing oj. For Rydberg states with n^l 
and low L one has (ci;L^/3)-Cl. In this case Eq. ^ can 
be simplified to the form: 



cr^ {nL ~> E,L±1) 
1 



9cn^ 



64/3r2(2/3) , 62/3r2(l/3) 2' 



,7/3 



,5/3 



(7) 



The combination of Eqs. (g]) and ^ yields: 



Wbbr 



2.22w"^/^-h2.63w^/^L^ 
duj 



QUj/{kT) _ I ■ 



(8) 



The expression in square brackets is a slowly varying 
function of uj. Taking into account that the main con- 
tribution to Wbbr is due to frequencies near the ion- 
ization threshold, one can replace uj by 1/(271^). After 
such replacement the integral in Eq. ([8]) can be calculated 
analytically, and the final result is: 



Wbbr 



kT 



2.80 2.09L'^ 



,7/3 



7,11/3 



In 



l_exp (-^) 



(9) 

Equation ([9]) gives the approximate direct BBR pho- 
toionization rate in atomic units for T measured in 
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BBR-induced photoionization rate of the sodium nS (a), nP (b), nD(c) and rubidium nS (d), nP (e) and nD 
n—5-80 calculated at the ambient temperatures T—77, 300 and 600 K. 



Kelvins. Alternatively, it can be rewritten to yield 
Wbbr in the units of s~^ for temperature T taken in 
Kelvins as follows: 



BBR 



= CtT 




(10) 



Here C l is an L-dependent scaling coefficient, which will 
be discussed later. By replacing n with the effective prin- 
cipal quantum number, Eq. (jlOp can be used for quick 
estimations of direct BBR ionization rates for various 
Rydberg atoms and states. The precise values, however, 
should be calculated numerically. 

The results of our numerical and analytical calcula- 
tions of the direct BBR-induced photoionization rates for 
sodium and rubidium nS, nP and riD states with 7i=5-80 
at the ambient temperatures T=77, 300 and 600 K are 
shown in Fig. 2. A good agreement of our numerical re- 
sults obtained using the Dyachkov and Pankratov model 
with the theoretical data obtained by Lehman [12] is ob- 
served. For the case of rubidium such comparison is not 
possible because no other published data are available, 
to the best of our knowledge. 

It is also interesting and instructive to compare the 
results obtained by numeric calculation with those ob- 
tained using the analytical formula (fTU]) . Figure 2(a) 
shows a remarkable disagreement between the numeri- 



cal and analytical results (assuming Cl = 1) for sodium 
nS states, which are known to possess a large quantum 
defect ((5s=1.348). At the same time. Figs. 2(b) and 2(c) 
show that in the case of Na nP and nD states, which have 
smaller quantum defects ((5p=0.855 and (5_d=0.015), the 
agreement between the numerical and analytical results 
is much better. This is not unexpected, since Eq. ^ 
is valid only for states with small quantum defects. For- 
mally, the disagreement for the non-hydrogenic nS states 
stems from peculiarities of the asymptotic behavior of 
Bcssel functions in Eq. ^ for states with i <C 1: the an- 
alytical expression of GDK model yields close photoion- 
ization cross-section values for nS, nP and nD states, 
while the accurate numerical calculations yield signifi- 
cantly smaller cross-sections for sodium nS states [see 
Fig. 2(a)]. At the same time, one can see from Fig. 2(a) 
that shapes of the analytical curves are quite similar to 
the numerical ones. Therefore, one can simply introduce 
a scaling coefficient in Eq. pTI)) in order to make it valid 
also for nS states. 

In order to illustrate that, in Figs. 2(d)-(f) we show the 
rescaled results of Eq. pU]) for the case of Rb Rydberg 
atoms nS, nP, and nD states, which all have large quan- 
tum defects {5l = 3.13, 2.66, and 1.34, respectively). A 
good agreement between the analytical and numerical re- 
sults was obtained when the rate obtained from Eq. ITUl) 
was scaled by a factor of C5=0.2, Cp=0.2 and Cd=0.4 
for the case of nS, nP, and nD states, respectively. 

Our precise numerical data on W b br are summarized 
in Tables (|J)-(|Vl| of Appendix. 



5 



B. BBR-induced mixing of Rydberg states 

BBR causes not only direct photoionization of the ini- 
tially populated levels. It also induces transitions be- 
tween neighboring Rydberg states, thus leading to a pop- 
ulation redistribution 0, For example, after laser 
excitation of the Na 16S state, the BBR-induced tran- 
sitions populate the neighboring n'P states [Fig. 1(a)]. 
The calculations show that these states have significantly 
higher direct photoionization rates W bbb. than the 16S 
state itself. Hence, BBR-induced population transfer to 
n'P state can noticeably affect the measured effective 
BBR ionization rate. The rates of spontaneous and BBR- 
induced transitions from the initial 16S and 16D states 
to a number of n'P states have been calculated by us in 
prj and are shown in Figs. 1(b) and 1(c). 

Importantly, absorption of BBR induces also transi- 
tions to higher Rydberg states, which are denoted as n" 
in Fig. 1(a). These states can be field ionized by the elec- 
tric field pulses usually applied in experiments in order 
to extract ions into channeltron or microchannel plate 
detectors. 



C. Field ionization of high Rydberg states 
populated by BBR 

Extraction electric field pulses, which are commonly 
used to extract ions from the ionization zone to the ion- 
ization detector, ionize Rydberg states with principal 
quantum numbers n exceeding some critical value ric- 
This critical value ric depends on the amplitude of the 
applied electric field and it can be found from the ap- 
proximate formula jlZ] 



E, 



3.2 ■ lO^n"" (V/cm) . 



(11) 



were Ec is the critical electric field for Uc- Hence, if a 
BBR mediated process populates a state with n' > Uc, 
this state will be ionized and thus will contribute to the 
detected ionization signal 0. 

We calculated the radial matrix elements 
R [nL n'L') of all dipole-allowcd transitions to 
other n'L' states with L' = (L ± 1) using the semi- 
classical formulae of [l^. The rate of a BBR-induced 
transition between the states nL and n'L' is expressed 
through the rate of spontaneous transitions given by the 
Einstein coefficient A {nL n'L'): 



A (nL n'L') - A^ll!l^ji2 ^ 



W (nL n'L') = A {nL n'L' 



gc^„„,/(feT) _ I 



, (12) 



where = l/(2n^) — l/(2n'^) is the transition fre- 

quency. Note that spontaneous transitions are possible 



only to lower levels, while BBR leads to transitions both 
upwards and downwards. 

The total rate Wsfi of BBR transitions to all Ryd- 
berg states with n' > n,, was calculated by summing the 
individual contributions of nL — » n'L' transitions given 
by Eq. 



Ws 



FI 



W{nL^n'L') 



(13) 



n'>nc L'—L±\ 



The values of Wsfi were numerically calculated for 
various amplitudes E of the electric field pulses. 

We also compared the numerical values with those ob- 
tained from the approximate analytical formulae, which 
has been derived with the bound-bound matrix elements 
of the GDK model: 



Wsfi 



l/(2n2)_i/(2„2) 



(14) 



The integration limits are chosen such that the inte- 
gral accounts for transitions to those Rydberg states, for 
which [l/(2n^) — l/(2n^)] < cu < ujnL (i.e., states above 
the field ionization threshold). Integration of Eq. 
gives another useful analytical formula that is similar to 
Eq. HOl): 



14423 10770L2 



Wsfi = ClT 



,7/3 



,11/3 



X In- 



1 — exp 



- In- 



f 1578 



7890 _ 157890 \ 



1 — exp I 



157890 ^ 



(15) 



where T is in Kelvins. 

The obtained numerical and analytical data on W sfi 
are presented in Fig. 3. For Na atoms, the calculations 
were made for the nS, nP and nD states with n=5-35 at 
the ambient temperature r=300 K [Figs. 3(a), 3(b) and 
3(c)]. The amplitudes E of the extracting electric field 
pulse are chosen as 100 and 200 V/cm (corresponding to 
nc = 42 and 36, respectively). These values are close 
to those used in our recent atomic-beam experiment on 
ionization of Na atoms llj. Alternatively, Figs. 3(d), 
3(e) and 3(f) show the calculated rates Wsfi for nS, nP 
and nD states of Rb for the lower E values (5 V/cm and 
10 V/cm corresponding to nc=91 and n^ =76, respec- 
tively), which are more adequate to experiments with 
Rydberg atoms in cold gases. 
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FIG. 3: The calculated BBR-induced SFI rate W sfi for (a) Na(nS), (b) Na(nP), and (c) Na(nD) states with n=5-35 at the 
electric field amplitudes of iJ=100 and 200 V/cm, and for the (d) Rb(nS), (e) Rb(nP), and (f) Rb(nD) states with n=5-80 at 
the electric field amplitudes of E=h and 10 V/cm. In all cases the ambient temperature is r=300 K. 



In order to achieve a better agreement between the an- 
alytical formula (fT5|l and the results of numerical calcula- 
tions, scaling coefficients Cl in Eq. pS]) were also intro- 
duced. The best agreement was obtained at Cs — 0.2, 
Cp = 0.15 and Cu = 0.3 for the Rb nS, nP and nD 
states, respectively. In the case of nP and nD states 
of Na data such scaling was not necessary, while in the 
case of nS states of Na the best agreement [not shown 
in Fig. 3(a)] was found at Cs — 0.25. The possibility to 
achieve a satisfactory agreement between the numerical 
and analytical data of Fig. 3 suggests that, using appro- 
priate scaling coefficients, the analytical formula (fTS)) is 
also suitable to quick estimates of BBR-induced SFI rates 
for various Rydberg atoms and states. 

Our precise numerical data on Wsfi are summarized 
in Tables (|J)-(|VT| of Appendix. 



by the laser pulse. Then atoms are allowed to interact 
with ambient BBR during the time interval (ti, ^2)- The 
second electric-field pulse extracts the ions, which have 
been produced by coUisional and BBR-induced ioniza- 
tion, to the ion detector. The atomic and molecular ions 
were distinguished using a time-of-flight method. At the 
ground-state density of 5 x 10^°cm~^ the BBR ionization 
is the main source of atomic ions [lT| and the contribution 
from coUisonal ionization of Penning type is negligible. 

Let us consider first the simplest case of laser excitation 
of a single nS state. The evolution of the number iV^+ 
of atomic ions produced via absorption of BBR photons 
by atoms in the initial nS state is given by 



d^A^ it) 
dt 



= WBBRNnS (t) , 



(16) 



D. Ionization of Rydberg states populated by BBR 

In this section we shall analyze the influence of time 
evolution of populations of Rydberg states upon interac- 
tion with ambient BBR photons. The typical timing dia- 
gram of laser excitation of Rydberg states and detection 
of ions is shown in Fig. 4. Such scheme was used in our 
recent experiment on collisional ionization of Na Rydberg 
atoms [llj. The first electric-field pulse is applied imme- 
diately after the laser excitation pulse in order to remove 
the atomic and molecular photoions produced 



where Nns{t) = N^sit = 0)exp{-t/T^ff) is the total 
number of Rydberg atoms remaining in the nS state as 
a function of time, and t^^j- is an effective lifetime of the 
nS state at a given ambient temperature. The registered 
photoions are produced during the time interval (ii,t2). 
The total number of ions can then be found by integrating 
Eq. HH) from ti to ij- 



Na+ =NnS it = 0)WBBRX 

[exp i^h/T^ff) - exp 



XT, 



nS 



(17) 
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FIG. 4: Timing diagram of signals: (a) laser excitation pulse; 
(b) electric field pulses for the ion extraction; (c) atomic A'^ 
and molecular ylj ion signals; (d) detector gates. 



This result can be rewritten by introducing an eflPective 
interaction time 11 ll: 



Na+ = Nns {t 



0) 



nS 



[exp {-h/r^ff) - exp {^h/r^f^)] . (18) 



Blackbody radiation also induces transitions to other 
Rydberg states n'P, as was discussed in Section II B. 
Evolution of populations of these states is described by 
the rate equation 



dN^'P jt) 
dt 



= [W {nS n'P) + A{nS^ n'P)] N^s {t) - 

-N^,P {t)lT%P, (19) 

where ^(nS n'P) and Ty(nS n'P) are the rates 
of population of the n'P state due to spontaneous tran- 
sitions to lower levels and BBR induced transitions up- 
wards and downwards from the initial nS state, and r"^^ 
is the effective lifetime of the n'P state. 

A combination of Eq. with Eqs. (HI]) and (fTT]) 
yields 
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FIG. 5: Calculated total BBR-induced ionization rates W^gn 
for (a) nS and (b) nD states of Na with n=8-35 for the ex- 
tracting electric field pulses of 100 V/cm and 200 V/cm at the 
ambient temperatures of T—77, 300 and 600 K. Experimental 
points are taken from Ref. |TH. 



states to highly excited states n"L" with n" > nc [see 
Fig. 5(a)], followed by ionization of these states by ex- 
tracting electric field pulses. This rate can be calculated 
using the same Eq. (f20|) . in which Wbbb. is replaced by 
WsFi and the summation is done over the states with 
n' < nc- 

Our precise numerical data on W^jfj^ and W™pj are 
summarized in Tables (II))- (jVI|) of Appendix. 



III. RESULTS OF CALCULATIONS AND 
COMPARISON WITH EXPERIMENT 



(c;,") -(-?//) 

xWBBRin'P) ( 1 • (20) 

The main contribution to the sum in Eq. (j20p is from 
n'P states with n' n ± 1 [see Fig. 1(b)]. The effec- 
tive BBR ionization rates for nP and nD states were 
determined in the same way as for nS states, taking into 
account the population transfer to both n'{L + 1) and 
n'{L — 1) states. 

The rate W™p^ describes the second-order process 
of BBR-induced transitions from the neighboring n' L' 



Calculated total BBR-induced ionization rates W^gj^ 
for Na nS, nP and nD states at the temperatures T = 
77, 300, 600 K and for the extracting electric field pulses 
of i?=100 v/cm and 200 V/cm, as well as their compari- 
son with our experimental data from Ref. pH , are shown 
in Fig. 5. Since the values of Wggj^ depend on the tim- 
ing of the experiment (see Fig. 4), the calculations were 
performed for <i ~ 0.3 /iS and t2 = 2.1 /is, which were 
used in our experiment JjJ. At the 300 K temperature 
a good agreement is observed for nD states with n=:8-20 
and for nS states with n=8-15. At the same time, for nS 
states with n > 15 the calculated ionization rates start 
to deviate from the experimental values; the measured 
values exceed the theoretical ones by a factor of 2.1 for 
n=20, and the shape of the experimental n-dependence 
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FIG. 6: The calculated total BBR-induced ionization rates Wbbr for (a) Na(nS), (b) Na(nP), (c) Na(nD), (d) Rb(nS), (e) 
Rb(nP), and (d) Rb(TiD) Rydberg states for the electric field amplitudes of 5 V/cm and 10 V/cm at the ambient temperatures 
of r=77 K, 300 K and 600 K. 



differs from the theoretical one. 

One possible explanation of such anomaly for nS states 
is related to their specific orbit that penetrates into the 
atomic core. The penetration causes a strong interac- 
tion between the Rydberg electron and the core, e.g., 
due to core polarization [16|]. This results in large quan- 
tum defect and a Cooper minimum in the photoionization 
cross-sections. This assumption is supported by the good 
agreement of theory and experiment for the hydrogen-like 
TiD states, which have a small quantum defect and almost 
non-penetrating orbits. 

Total BBR-induced ionization rates were also calcu- 
lated for Na and Rb nS, nP and nD states in a broader 
range of n and for lower amplitudes of the electric field 
pulses (5 V/cm and 10 V/cm). Such fields are more rel- 
evant to the experiments with cold Rydberg atoms, e.g., 
to ultraco ld p lasma formation from cold Rydberg atoms 
clouds [§, [ij , since these experiments explore the atoms 
excited to relatively high Rydberg states that are ion- 
ized in weaker electric fields (35 V/cm for n 50 and 
10 V/cm for n ~ 80). The calculated ionization rates 
are presented in Fig. 6. It can be seen that SFI and 
BBR-induced level-mixing processes alter the shapes of 
n-dependences of total ionization rates Wggj^, which ef- 
fect is more pronounced at lower temperatures and larger 
n. As Fig. 6 has a logarithmic scale, the visible effect of 
additional processes looks small. However, from the ta- 
bles of Appendix it can be seen that, e.g., for the Rb 65S 
state at T=300 K the total BBR ionization rate is twice 
larger than the direct photoionization rate (see Table 5). 



Precise results of our numerical calculations for a five- 
step set of principal quantum numbers n are summarized 
in Tables (|I)l- (jVip of Appendix. The values of ionization 
rates at the other n values can be obtained either by 
extrapolation or upon request from the authors of this 
paper. 

IV. CONCLUSION 

We have calculated the total BBR-induced ionization 
rates of Na and Rb nS, nP and nD Rydberg states for 
principal quantum numbers n=8-65 at the ambient tem- 
peratures of 77, 300 and 600 K. Our calculations take 
into account the effect of BBR-induced mixing of Ryd- 
berg states and their field ionization by extracting elec- 
tric field pulses. Useful analytical formulae have been 
derived, which allow for quick estimation of ionization 
rates and their dependences on the principal quantum 
number n. The numerical results are in a good agree- 
ment with our recent experiment data on Na nS and nD 
states, except for nS states with n > 15, which is most 
probably associated with the Cooper minimum in the 
photoionization cross-section. 

The obtained results show that BBR-induced redis- 
tribution of population over Rydberg states and their 
field ionization by extracting electric fields affect both 
the magnitudes of total ionization rates and shapes of 
their dependencies on the principal quantum number. 
This suggests that these processes are important and can- 



not be ignored in the calculations and measurements of 
BBR ionization rates. Equations P^ - (|20p . as well as the 
analytical formulae Q and (fT5|l . can be used to calcu- 
late total ionization rates W^gj^ under particular exper- 
imental conditions. The obtained numerical results may 
be helpful to the analysis of ionization signals measured 
in experiments on coUisional ionization and spontaneous 
formation of ultracold plasma, since BBR-induced ion- 
ization is the main channel of delivering atomic ions. At 
the same time, as we have revealed that theoretical data 
for Na rtS-states noticeably disagree with experiment at 
15 < n < 20, new experimental data for Na and Rb in a 
broader range of principal quantum numbers would be of 
interest for the further improvement of theory, especially 
for the non-hydrogen-like states. 
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TABLE I: Calculated BBR-induced ionization rates (s ^) for Na at T=77 K 



nS 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 




w^M 

lOOV/cm 


200V/cm 


W'b''br 
lOOV/cm 


Wj3°BR 

200V/cm 


10 


1.56 X 10"^ 


3.44 X 10"^ 


7.71 X 10"^ 


4.92 X 10"^* 


3.55 X 10"'' 


7.03 X 10"'' 


2.12 X lO"'^ 


2.2 X 10"'' 


15 


0.009 


0.019 


0.042 


0.002 


0.003 


0.006 


0.033 


0.059 


20 


0.805 


1.655 


3.83 


0.059 


0.118 


0.240 


2.64 


4.94 


25 


4.04 


10.23 


25.29 


0.185 


0.420 


0.891 


14.9 


30.4 


30 


8.45 


27.12 


81.09 


0.257 


0.706 


1.689 


36.5 


91.5 




Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cm 


Tirmix 
BBR 


TjTjnix 
VVsFI 

5V/cm 


Til/mix 

lOV/cm 


VV BBR 

5V/cm 


VV BBR 

lOV/cm 


35 


11.8 


2.47 


5.20 


0.267 


0.0488 


0.107 


14.6 


17.4 


45 


14.6 


4.12 


9.09 


0.195 


0.0466 


0.106 


19.0 


24.0 


55 


14.1 


5.76 


13.9 


0.124 


0.0412 


0.102 


20.1 


28.3 


65 


12. () 


8.01 


2:!.() 


0.0773 


0.0377 


0.111 


20.7 


:!(). 1 


nP 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


BBR 


VV spi 

lOOV/cm 


VVgpi 

200V/cm 


Wh%R 
lOOV/cm 


Wh%R 
200V/cm 


10 


7.09 X 10"® 


3.49 X 10"'' 


6.91 X lO"'^ 


4.75 X 10"® 


1.51 X 10"® 


2.99 X 10"® 


7.5 X 10"® 


7.85 X 10"® 


15 


0.0996 


0.195 


0.389 


0.00111 


0.00200 


0.00409 


0.298 


0.493 


20 


5.22 


10.2 


20.5 


0.0217 


0.0455 


0.0977 


15.4 


25.8 


25 


22.1 


50.8 


107 


0.0545 


0.132 


0.308 


73.1 


130 


30 


12.1 


120 


290 


().0()()2 


0.197 


0.551 


Ki:! 






Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cm 


vvrbr 


■tirmix 
VV spi 

5V/cm 


lOV/cm 


Wh%R 
5V/cm 


i^rtot 
Vy BBR 

lOV/cm 


35 


57.0 


10.5 


22.9 


0.0726 


0.0144 


0.0318 


67.6 


80.0 


45 


(i(i.2 


1().2 


:!7.1 


0.0511 


().01:i5 


0.():ii:i 


82.5 


10:! 


55 


61.9 


21.6 


53.7 


0.0320 


0.0117 


0.0297 


83.5 


116 


65 


53.8 


28.6 


86.5 


0.0186 


0.00997 


0.0304 


82.4 


140 


nD 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


BBR 


Tirmix 
VVSFI 

lOOV/cm 


Wspi 

200V/cm 


W*B°U 

lOOV/cm 


Wh°*BR 

200V/cm 


10 


1.11 X 10"" 


1.53 X 10"^ 


2.98 X lO"'^ 


7.61 X lO"'^ 


4.96 X lO"'^ 


9.82 X 10"^ 


1.28 X 10"" 


1.43 X 10""* 


15 


0.251 


0.483 


0.947 


0.00270 


0.00472 


0.00953 


0.742 


1.21 


20 


G.59 


12.7 


25.2 


0.0371 


().07()2 


0.158 


19.1 


32.0 


25 


22.6 


50.0 


104 


0.0926 


0.215 


0.467 


72.9 


127 


30 


39.0 


107 


256 


0.120 


0.331 


0.812 


146 


296 




Wbbr 
direct 


WsFI 

.")\'/'cm 


10\'/cm 


1' BBR 


V\ SFI 

5\'/'cm 


Tjr7nix 
VV SFI 

l()\'/(:m 


BBR 
5\'/(:m 


Tirtot 

BBR 

l()\'/cm 


35 


49.3 


8.88 


19.4 


0.131 


0.0239 


0.0527 


58.3 


68.8 


45 


54.6 


13.0 


29.6 


0.0995 


0.0230 


0.0528 


67.7 


84.3 


55 


49.8 


16.7 


41.5 


0.0663 


0.0204 


0.0503 


66.7 


91.5 


65 


42.7 


21.8 


65.9 


0.0424 


0.0180 


0.0519 


64.6 


109 
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TABLE II: Calculated BBR-iuchiccd ionization rates (s^^) for Na at r=300 K 



nS 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 




Ws'fJ 
lOOV/cm 


200V/cm 


Wbbr 
lOOV/cm 


Wbbr 
200V/cm 


10 


7.50 


1.33 


2.34 


13.6 


3.02 


4.88 


25.5 


28.3 


15 


114 


39.4 


70.2 


78.6 


24.2 


39.6 


256 


302 


20 


223 


99.2 


185 


77.6 


29.2 


49.2 


429 


534 


25 


254 


155 


314 


53.6 


26.2 


46.9 


489 


669 


30 


242 


218 


537 


34.2 


23.2 


46.9 


518 


861 




Wbbr 
direct 


WSFI 

5V/cm 


WsFI 

lOV/cm 


Wbbr. 


WSFI 

5V/cm 


WsFI 

lOV/cm 


W'b%r 
5V/cm 


Wbbr 
lOV/cm 


35 


217 


19.3 


39.5 


22.6 


1.55 


3.29 


261 


283 


45 


166 


22.8 


49.1 


10.2 


1.04 


2.30 


200 


227 


55 


126 


27.3 


64.2 


5.15 


0.780 


1.87 


159 


197 


65 


97.4 


35.0 


100 


2.80 


0.652 


1.86 


136 


202 


nP 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


Wbbr 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


W'B°^a 
lOOV/cm 


Txrtot 

Wbbr 
200V/cm 


10 


79.92 


15.99 


25.84 


2.426 


0.8949 


1.456 


99.23 


109.6 


15 


778.2 


240.9 


393.6 


17.82 


5.476 


9.174 


1042 


1199 


20 


1295 


501.8 


845.4 


19.05 


6.384 


11.24 


1823 


2171 


25 


1374 


709.0 


1272 


13.50 


5.610 


10.76 


2102 


2670 


30 


125:! 


9:52. 


1919 


8.111 


1. ()()() 


10.56 


2198 


:5191 




Wbbr 
direct 


WSFI 

5V/cm 


WSFI 

lOV/cm 


\Jrrnix 

Wbbr 


WsFI 

5V/cm 


WSFI 

lOV/cm 


W'b%r 
5V/cm 


lOV/crn 


35 


1090 


79.97 


169.6 


6.145 


0.3428 


0.7388 


1177 


1267 


45 


795.2 


89.2(S 


198.3 


2.830 


0.2290 


0.5203 


887.5 


996.9 


55 


586.0 


102.1 


247.0 


1.419 


0.1679 


0.4172 


689.7 


834.8 


65 


441.8 


124.9 


366.0 


0.7254 


0.1267 


0.3768 


567.5 


808.9 


nD 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


Wbbr 


lOOV/cm 


WsFI 

200V/cm 


Wbbr 
lOOV/cm 


Tirtot 

Wbbr 
200V/cm 


10 


157.5 


37.01 


59.22 


13.71 


2.596 


4.237 


210.8 


234.7 


15 


890.1 


268.2 


434.4 


39.52 


8.003 


13.30 


1206 


1377 


20 


1253 


467.7 


782.0 


.36.76 


8.176 


14.05 


1765 


2086 


25 


1211 


013. 5 


1091 


2().17 


6.907 


12.61 


1891 


2:!77 


30 


1099 


774.0 


1592 


17.38 


5.679 


11.70 


1896 


2720 




Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cin 


Wbbr. 


5V/ciii 


lOV/cin 


wh%R 

5V/crn 


lOV/cin 


35 


936.2 


()."). 00 


i:i7.7 


12. .31 


0.1207 


0.8980 


1011 


1087 


45 


666.5 


70.18 


155.6 


6.040 


0.2783 


0.6197 


743.0 


828.8 


55 


483.1 


78.52 


189.6 


3.240 


0.2025 


0.4856 


565.1 


676.4 


65 


360.0 


94.60 


277.3 


1.849 


0.1557 


0.4332 


456.6 


639.6 
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TABLE III: Calculated BBR-induced ionization rates (s"^) for Na at r=600 K 



nS 


n 


Wbbr 
direct 


WsFI 

lOOV/crn 


WsFI 

200V/cm 




w^M 

lOOV/cm 


w^M 

200V/cm 


lOOV/cm 


W'bsh 
200V/cm 


10 


284.4 


41.90 


71.08 


1396 


187.6 


294.9 


1910 


2046 


15 


910.1 


185.6 


321.5 


1396 


237.4 


378.9 


2729 


3007 


20 


991.5 


290.5 


526.8 


770.7 


175.8 


289.9 


2229 


2579 


25 


874.7 


378.8 


748.8 


411.1 


130.9 


229.0 


1796 


2264 


30 


728.5 


489.6 


1172 


229.0 


105.5 


208.4 


1553 


2338 




Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cm 


TT/mix 
BBR 


WsFI 

5V/cm 


Tjrmix 
WsFI 

lOV/cm 


w%°U 

5V/cm 


yy BBR 

lOV/cm 


35 


601 


42.9 


87.6 


139 


7.04 


14.8 


790 


843 


45 


416 


48.3 


104 


57.1 


4.46 


9.82 


526 


586 


55 


299 


56.5 


132 


27.5 


3.25 


7.76 


386 


467 


65 


222 


71.3 


2();i 


1 l.G 


2.07 


7.57 


311 


118 


nP 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


BBR 


■i^rmix 
"SFI 

lOOV/cm 


200V/cm 


Wh%R 
lOOV/cm 


wh°hR 

200V/cm 


10 


2580 


353.0 


554.1 


190.1 


31.47 


49.95 


3155 


3374 


15 


6054 


1054 


1681 


279.2 


45.01 


73.60 


7433 


8088 


20 


5900 


1433 


2364 


173.8 


34.97 


60.14 


7541 


8497 


25 


4884 


1715 


3015 


96.59 


26.01 


48.63 


6722 


8044 


30 


3898 


2077 


4183 


52.98 


19.78 


43.40 


6048 


8178 




Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cm 


y^ BBR 


Wmix 
W spi 

5V/cm 


ijTmix 
W SFI 

lOV/cm 


yy BBR 

5V/cm 


y^ BBR 

lOV/cm 


35 


3115 


177.6 


375.1 


35.25 


1.458 


3.135 


3330 


3529 


45 


2049 


188.9 


417.8 


14.72 


0.9238 


2.094 


22.54 


2484 


55 


i 12() 


210.8 


508.0 


7.008 


().()582 


l.():i2 


Kill 


1912 


65 


1029 


254.5 


742.4 


3.470 


0.4856 


1.439 


1288 


1777 


35 


3115 


177.6 


375.1 


35.25 


1.458 


3.135 


3330 


3529 


nD 


n 


Wbbr 
direct 


WsFI 

lOOV/cm 


WsFI 

200V/cm 


TJ/mix 
" BBR 


VVSFI 

lOOV/cm 


Ti/mia; 
W spi 

200V/cm 


yy BBR 

lOOV/cm 


" BBR 

200V/cm 


10 


3679 


506.3 


787.9 


750.7 


68.35 


108.5 


5004 


5326 


15 


6226 


1052 


1666 


631.2 


58.47 


94.75 


7968 


8618 


20 


5547 


1286 


2106 


.368.8 


40.12 


67.39 


7242 


8090 


25 


1105 


1 158 


2549 


211.4 


28.43 


50.80 


0103 


7210 


30 


3434 


1709 


3438 


124.9 


21.11 


42.54 


5288 


7039 




Wbbr 
direct 


WSFI 

5V/ciii 


WSFI 

lOV/ciii 


WTbr 


•iirmix 
WsFI 

5V/cni 


XT rmix 

lOV/cm 


wh°hn 

5V/cm 


Wbbr 
lOV/crn 


35 


2702 


113.5 


:i()2.9 


81.()() 


1.591 


:!.:i9() 


2929 


301)0 


45 


1733 


148.1 


327.0 


36.94 


0.9921 


2.202 


1919 


2099 


55 


1185 


162.0 


389.4 


19.08 


0.6971 


1.667 


1367 


1595 


65 


844.1 


192.7 


562.0 


10.72 


0.5188 


1.438 


1048 


1418 
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TABLE IV: Calculated BBR-induced ionization rates (s"^) for Rb at T=77 K 



nS 


n 


Wbbr 
direct 


WsFI 

5V/cm 


WsFI 

lOV/cm 


TJ/mix 
BBR 


T I rmix 
SFI 

5V/cm 


WsFI 

lOV/cm 


W'b°'br 
5V/cm 


W'b^'br 
lOV/cm 


10 


3.34 X 10"^^ 


4.54 X 10"^^ 


9.46 X 10"" 


3.36 X 10"" 


4.36 X 10"" 


9.21 X 10"" 


4.18 X 10"" 


4.72 X 10"" 


20 


0.2826 


0.04174 


0.08730 


0.005807 


7.58 X 10^* 


0.001622 


0.3309 


0.3773 


30 


6.541 


1.108 


2.343 


0.04932 


0.006986 


0.01519 


7.706 


8.949 


40 


12.61 


2.740 


5.967 


0.04750 


0.008213 


0.01840 


15.41 


18.64 


50 


13.91 


4.201 


9.678 


0.03120 


0.007135 


0.01687 


18.15 


23.64 


60 


12.97 


5.813 


14.96 


0.01917 


0.006156 


0.01615 


18.81 


27.97 


65 


12.21 


6.863 


19.53 


0.01502 


0.005856 


0.01688 


19.10 


31.78 


nP 


10 


2.96 X 10"^° 


3.84 X 10"" 


8.1 X 10"" 


2.51 X 10"" 


3.39 X 10"" 


7.1 X 10"" 


3.35 X 10"" 


3.77 X 10"" 


20 


0.489 


0.0636 


0.136 


0.00443 


6.31 X 10"" 


0.00132 


0.558 


0.631 


30 


7.21 


1.02 


2.22 


0.0412 


0.00656 


0.0139 


8.28 


9.49 


40 


11.8 


2.08 


4.65 


0.0447 


0.00887 


0.0193 


14.0 


16.5 


50 


11.9 


2.83 


6.70 


0.0321 


0.00854 


0.0195 


14.8 


18.7 


60 


10.5 


3.60 


9.49 


0.0212 


0.00797 


0.0201 


14.1 


20.0 


65 


9.65 


4.09 


11.9 


0.0171 


0.00785 


0.0217 


13.8 


21.6 


nD 


10 


5.11 X 10"'^ 


6.76 X 10"® 


1.42 X 10"'' 


4.23 X 10"^ 


5.57 X 10"" 


1.18 X 10"® 


5.83 X lO"'^ 


6.58 X 10"^ 


20 


1.66 


0.227 


0.475 


0.0149 


0.00207 


0.00442 


1.90 


2.15 


30 


15.7 


2.46 


5.19 


0.0710 


0.0110 


0.0237 


18.2 


21.0 


40 


24.7 


4.99 


10.8 


0.0662 


0.0127 


0.0280 


29.7 


35.6 


50 


25.4 


7.15 


16.5 


0.0450 


0.0113 


0.0264 


32.6 


41.9 


60 


22.9 


9.60 


24.9 


0.0287 


0.0101 


0.0260 


32.5 


47.8 


65 


21.3 


11.3 


32.5 


0.0229 


0.00969 


0.0277 


32.6 


53.9 
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TABLE V: Calculated BBR-induced ionization rates (s"^) for Rb at r=300 K 



nS 


71 


vVbBR 

direct 


vVSFI 

5V/cm 


vVSFI 

lOV/cm 


>^BBR 


Wspj 

5V/cm 


Wspj 
lOV/cm 


TT/to£ 

Wbbr 
5V/cm 


Ti/to£ 

Wbbr 
lOV/cm 


10 


0.604 


0.0195 


0.0390 


0.408 


0.0113 


0.0228 


1.04 


1.07 


20 


183 


7.58 


15.3 


20.8 


0.633 


1.30 


212 


220 


30 


235 


14.8 


30.3 


9.25 


0.390 


0.819 


260 


276 


40 


189 


18.5 


39.2 


3.70 


0.227 


0.492 


211 


232 


50 


143 


21.9 


49.1 


1.67 


0.150 


0.344 


167 


194 


60 


109 


26.7 


67.0 


0.842 


0.113 


0.288 


137 


177 


65 


96.1 


30.3 


83.9 


0.619 


0.103 


0.287 


127 


181 


nP 


10 


3.69 


0.102 


0.206 


0.0598 


0.00432 


0.00873 


3.85 


3.96 


20 


:!()() 


9.25 


19.0 


8.11 


0.:ill 


0.()2() 


:!17 


:!27 


30 


294 


13.0 


27.3 


4.87 


0.255 


0.520 


312 


327 


40 


207 


13.8 


30.0 


2.20 


0.173 


0.363 


223 


240 


50 


144 


14.7 


33.8 


1.07 


0.126 


0.278 


160 


179 


60 


103 


16.5 


42.3 


0.568 


0.103 


0.246 


120 


146 


65 


87.(S 


1«.0 


51.0 


I). 12() 


().09()5 


0.251 


lOG 


110 


nD 


10 


27.8 


0.782 


1.58 


2.73 


0.0532 


0.108 


31.4 


32.2 


20 


540 


18.7 


37.6 


21.4 


0.626 


1.28 


581 


600 


30 


517 


27.3 


55.9 


10.4 


0.439 


0.908 


555 


584 


40 


381 


31.6 


66.9 


4.60 


0.282 


0.600 


418 


453 


50 


277 


36.2 


81.2 


2.25 


0.199 


0.444 


315 


360 


60 


205 


43.5 


110 


1.20 


0.157 


0.385 


250 


316 


65 


179 


49.1 


138 


0.904 


0.146 


0.391 


229 


318 



TABLE VI: Calculated BBR-induced ionization rates (s^^) for Rb at r=600 K 



nS 


n 


Wbbr 
direct 


WsFI 

5V/crn 


WsFI 

lOV/crn 


Wbbr 


Wgpj 

5V/crri 


Wgpj 

lOV/cni 


Wbbr 
5V/cm 


w bbr 
lOV/cm 


10 


81.20 


1.429 


2.838 


177.4 


2.395 


4.806 


262.4 


266.2 


20 


941.6 


26.22 


52.51 


259.5 


4.569 


9.340 


1232 


1263 


30 


734.1 


35.55 


72.58 


68.99 


1.925 


4.022 


840.6 


879.7 


40 


500.8 


40.35 


85.19 


23.30 


1.009 


2.181 


565.4 


611.4 


50 


350.5 


45.92 


102.7 


9.740 


0.6394 


1.458 


406.8 


464.4 


60 


255.3 


54.93 


137.1 


4.738 


0.4709 


1.188 


315.4 


398.3 


65 


220.8 


61.92 


170.6 


3.445 


0.4253 


1.176 


286.6 


396.0 


nP 


10 


373.9 


5.011 


10.06 


20.19 


0.6466 


1.299 


399.7 


405.5 


20 


1690 


30.95 


63.33 


100.7 


2.073 


4.156 


1824 


1858 


30 


1009 


30.95 


64.80 


35.35 


1.170 


2.375 


1077 


1112 


40 


592.4 


30.00 


64.98 


13.44 


0.7148 


1.488 


636.5 


672.3 


50 


374.1 


30.71 


70.43 


6.029 


0.4986 


1.086 


411.3 


451.6 


60 


252.3 


33.85 


86.50 


3.076 


0.3940 


0.9362 


289.6 


342.8 


65 


211.3 


36.78 


103.7 


2.279 


0.3675 


0.9453 


250.7 


318.2 


nD 


10 


12()3 


17. 15 


:! 1.98 


279.7 


2.919 


5.900 


15():! 


158 1 


20 


2696 


57.56 


115.5 


220.0 


3.948 


8.007 


2978 


3040 


30 


1681 


63.97 


130.5 


71.48 


2.059 


4.235 


1819 


1887 


40 


1051 


68.29 


144.0 


27.59 


1.211 


2.556 


1148 


1225 


50 


699.6 


75.56 


168.9 


12.59 


0.8201 


1.814 


788.6 


883.0 


60 


192.7 


89.13 


223.8 


G.510 


0.0320 


1.534 


589.0 


724.5 


65 


420.5 


100.1 


280.1 


4.847 


0.5827 


1.544 


526.1 


706.9 
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